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Abstract 

The so-called parity-time-reversal- {VT-) symmetric quantum mechanics (PTQM) has developed 
into a noteworthy area of research. However, to date most known studies of PTQM focused on 
the spectral properties of non-Hermitian Hamiltonian operators. In this work, we propose an 
axiom in PTQM in order to study general time-dependent problems in PTQM, e.g., those with 
a time-dependent PT-symmetric Hamiltonian and with a time-dependent metric. We illuminate 
our proposal by examining a proper mapping from a time-dependent Schrodinger-like equation of 
motion for PTQM to the familiar timc-dcpcndcnt Schrodingcr equation in conventional quantum 
mechanics. The rich structure of the proper mapping hints that time-dependent PTQM can be a 
fruitful extension of conventional quantum mechanics. Under our proposed framework, we further 
study in detail the Berry phase generation in a class of T^T-symmetric two-level systems. It is 
found that a closed adiabatic path in PTQM is often associated with an open adiabatic path in a 
properly mapped problem in conventional quantum mechanics. In one interesting case we further 
interpret the Berry phase as the flux of a continuously tunable fictitious magnetic monopole, 
thus highlighting the difference between PTQM and conventional quantum mechanics despite the 
existence of a proper mapping between them. 

PACS numbers: 03.65.-w, 03.65.Vf, ll.30.Er 



I. INTRODUCTION 



As a fruitful extension of conventional quantum mechanics (QM), the so-called parity- 
time-reversal- iVT-) symmetric quantum mechanics (PTQM) has developed into a note- 
worthy area of research [HE]. So far, the main emphasis of PTQM studies has been placed, 
both theoretically and experimentally, on the spectral properties of non-Hermitian Hamilto- 
nians. Due to a nonconventional and well-behaved inner-product structure, the spectrum of 
a complex but PT-symmetric non-Hermitian Hamiltonian is still real if the VT symmetry 
is not spontaneously broken. This fascinating result hints the potential of introducing a 
non-trivial inner-product structure to QM. 

In conventional QM, the spectral problem for a Hamiltonian or energy operator is tackled 
by the stationary Schrodinger equation. The general equation of motion is given by a sec- 
ond equation, i.e., the time-dependent Schrodinger equation that describes how a quantum 
system evolves with time. There are opposite views about which of the two Schrodinger 
equations is more vital in conventional QM. For example, by considering whether an equa- 
tion is invariant upon a unitary transformation, the stationary Schrodinger equation seems 
more appealing. However, the stationary Schrodinger equation can be also regarded as a 
practical tool to solve the time-dependent Schrodinger equation, with the latter being a 
postulate in QM. Indeed, Dirac directly put down the time-dependent Schrodinger equation 
with an arbitrary time-dependent Hamiltonian, by reasoning that the Hamiltonian or the 
energy operator of a system should be a generating operator of time displacements [5]. 

Given the importance of the time-dependent Schrodinger equation in conventional QM, it 
is necessary to seek a general equation of motion in PTQM. The sought equation represents 
a time-dependent Schrodinger-like equation with a nonconventional inner product structure 
in Hilbert space. Such an equation of motion should be as general as possible, so as to cover 
cases with (i) time-dependent non-Hermitian Hamiltonians and (ii) time-dependent inner 
product structure. Efforts devoted to time evolution in PTQM will supplement ongoing 
PTQM activities on the spectral side. Time evolution in PTQM is also a necessary subject 
in order to treat conventional QM as a special case of PTQM, 

In our previous work, it was noted that for PTQM, there is ambiguity when constructing 
a time-dependent Schrodinger-like equation yielding unitary evolution [4J. In retrospect, 
this must be the case because even in conventional QM, the unitarity requirement does 
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not suffice to determine the form of a time-dependent Schrodinger equation. In this work, 
an axiom is proposed to remove the ambiguity. Our task is achieved by ffist noticing a 
dual role of Hamiltonian operators in conventional QM and then make a close analogy 
in PTQM. Further, we illuminate our axiom by examining the mapping of equation of 
motion from PTQM to that in conventional QM. We are able to identify a special type 
of mapping, termed as "proper mapping" below, from which we gain important insights 
into time-dependent PTQM. We use the complex harmonic oscillator problem as a simple 
working example to illustrate a proper mapping. 

To further elaborate our theoretical proposal, we revisit the Berry-phase problem for a 
class of PT-symmetric two-level systems Detailed Berry-phase calculations are pre- 
sented. We find that the results can be interpreted as the existence of a fictitious nontrivial 
Dirac string and a fictitious continuously tunable monopole. This theoretical finding is 
unexpected based on our intuition from conventional QM. We explain it by investigating a 
proper mapping between conventional QM and PTQM. In particular, it is found that a cyclic 
parameter path in PTQM is not necessarily a cyclic one with the same period for a prop- 
erly mapped problem in conventional QM. It is our belief that these specific results should 
stimulate further general interests in time-dependent PTQM, with our proposed axiom as a 
promising starting point. 

The organization of this paper is as follows. In Sec. |Tl], we first make remarks on the 
time-dependent Schrodinger equation in conventional QM and then propose a Schrodinger- 
like equation of motion for PTQM. We further define a proper mapping between these two 



theories. In Sec. |III| we adopt a simple example to illustrate our time-dependent PTQM. 
Detailed results of Berry phase calculations for a class of PT-symmetric two-level systems 



are presented and discussed in depth in Sec. [IV} Section |V] concludes this work. Some helpful 
derivations are also presented in Appendix A. 
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II. SCHRODINGER EQUATION OF MOTION: FROM CONVENTIONAL QM 
TO PTQM 



A. Remarks on Hamiltonian operators in conventional QM 

In conventional quantum mechanics, the Hamiltonian operator h determines the spectrum 
of a system through the stationary Schrodinger equation, 

/i|</)„) = (1) 

where £"„ are the eigenvalues and |0„) represent the eigenkets. Even if h is time-dependent, 
the above eigenvalue equation is still useful as it yields an instantaneous spectrum for the 
system (the instantaneous spectrum also gives the possible measurement results if the energy 
value is measured). Note that upon a unitary transformation U, the above equation is 
invariant with h ^ h' = UhW and — ^ \4>'n) = U\(f)n)- 

The time evolution of the above system is governed by the following time- dependent 
Schrodinger equation, 

ih^^mt)) = Ht)\m), (2) 

where we have used h{t) to stress that the Hamiltonian operator h can be explicitly time- 
dependent. Dirac justified this equation of motion by noting that, in this form the energy 
operator becomes the generator of time displacements, which is then analogous to the con- 
jugate pair of position and momentum variables [3]. However, such an equation of motion, 
which is a postulate in conventional QM, is only invariant under a time-independent unitary 
transformation. It is not invariant under arbitrary time- dependent unitary transformations. 
Indeed, if we apply a time-dependent unitary transformation, h{t) — )■ h'{t) = U{t)h{t)W{t) 
and ^ |$'(t)) = U{t)\^{t)), then Eq. Q becomes 

^\t)), (3) 

where an overhead dot represents the time derivative. As a result, some appropriate represen- 
tation must be implicitly adopted first when writing down the time-dependent Schrodinger 
equation Eq. (|2]). For our later use, we stress the evident observation here: the Hamiltonian 
or the energy operator in conventional QM plays a dual role: the operator that determines 
the energy spectrum of the system is also taken as the generator of time displacements. Later 



ih-\<l>'{t)) = h'{t)-ihU{t)U\t) 
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we shall exploit this dual role as a basis to define a "proper" mapping between conventional 
QM and PTQM. 



B. Time-dependent PT-symmetric quantum mechanics 

As is well known now, the Hilbert space in PTQM is associated with a metric operator W 
that defines an inner product We consider PTQM in its most general form by studying 
time-dependent Hamiltonian operators H{t) and also time-dependent metric operator W{t) 
associated with H{t). At any instant t under investigation, H{t) is assumed to have an 
extended VT symmetry characterized by 

W{t)H{t) = H{tyW{t). (4) 

The inner product between bra and ket states becomes (-IW^I-), where a bra state is defined 
as the Dirac conjugate of a ket state, (-1 = \-y . Note that we do not need to introduce a 
bi-orthonormal basis as commonly used in the study of open quantum systems with non- 
Hermitian Hamiltonians. The reason is that our PT-symmetric Hamiltonian is already 
assumed to be self-adjoint (under the metric W) when VT symmetry is unbroken. For a 
self-adjoint operator H{t), all its eigenvalues are real and its instantaneous eigenstates in 

H{t)\^4t)) = En{t)\Mt)) (5) 
form a set of orthonormal basis with respect to the metric W{t), 

{^m{t)\W{t)\ijrr{t))=S^n. (6) 

To construct an equation of motion, we impose the following unitary condition 

^^{^,{t)\wm2{t)) = o, (7) 

where |\E'i(t)) and |\E'2(i)) are two time evolving states from two different initial conditions 
|\l/i(0)) and |\l/2(0)). As already shown in our previous work [1], to satisfy the above unitary 
condition the equation of motion must take the following form 
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where the time-displacement generator A satisfies (from now on, we sometimes suppress the 
time argument) 

ihW = A^W - WA. (9) 

Interestingly, A must not equal H when the metric operator is changing in time [because 
WH — H^W = according to Eq. Q]. This hints that general time-dependent problems in 
PTQM go far beyond the spectral properties of a PT-symmetric Hamiltonian operator. 

Just like the fact that a generic matrix can be written as the sum of a Hermitian matrix 
and an anti-Hermitian matrix, we may also divide A into two parts, 

A = H + A, (10) 

with 

WH = H^W, (11) 

WA = -A^W. (12) 

Clearly, H is "PT-symmetric" under the same metric W as H. We call the second part A 
"PT-anti-symmetric" with respect to the metric operator W. In Appendix A, we show that 
for any given A and W, this partition is unique. 

Using the partition in Eq. (10), we rewrite Eq. ^ into 



ihW = A^W -WA = ~2WA. (13) 

This directly gives 

A = - lihW-^W. (14) 

Therefore, for a known W{t), the "PT-anti-symmetric" part of a general time-displacement 
generator A can be explicitly worked out. Summarizing, for a time-dependent PT-symmetric 
Hamiltonian H{t) which is diagonalizable with real eigenvalues (e.g., within a certain time 
window before PT-symmetry breaking), one may find a time-dependent metric W{t), and 
the generator of time displacements 

A = H- \ihW-^W (15) 

that yields unitary evolution, where H is always PT-symmetric with respect to W{t). 

It is important to stress that the above unitary condition does not lead to any conclusion 
about the relation between H and H (other than that they share the same metric W). This 



is not a surprise. Indeed, in convention QM, any Hermitian time-displacement generator 
leads to a unitary evolution. Taking the Hamiltonian operator as the generator of time 
displacements in conventional QM is a postulate, not a result by deduction. With this 
recognition, we proceed to propose an axiom for time-dependent PTQM, i.e., 

H = H, (16) 

which leads to 

A = H- \ihW-^W. (17) 

That is, the VT -symmetric part of the generator of time displacements in PTQM should he 
the same Hamiltonian operator that gives the (instantaneous) energy spectrum of the same 
system. We finally have the following time-dependent Schrodinger-like equation for PTQM, 

m-^l^) = {h- ^^W-^W^ I^). (18) 

This axiom lifts the ambiguity issue in PTQM, first raised in our previous work It can 
be seen that conventional QM can now be regarded as a special case of PTQM with the 
metric operator chosen to be unity at all times. 



C. Mapping between conventional QM and time-dependent PTQM 

Since both conventional QM and our time-dependent PTQM generate unitary evolutions, 
both yield real instantaneous spectra, one may wonder whether there is a simple mapping 
between these two frameworks. If there is such a mapping, then why do we still need 
time-dependent PTQM? In this subsection we shed light on these questions. 

As a metric operator, is a positive-definite Hermitian operator and hence it can be 
written as 

W = ri^f], (19) 
where t] is reversible. Note that the solution to r] is not unique and as a matter of fact, 

W = r]^r] = {Ut])^ {Ut]) with = (20) 

So any known solution multiplied by an arbitrary unitary operator (unitary in the sense 
of conventional QM) from the left is still a solution to rj. One may now use a known r) 
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to map a PT-symmetric operator H to a Hermitian operator h by the following similarity 
transformation, 

h = r]Hr]-\ (21) 



with h = being a direct consequence of the PTsymmetry of H [see Eq. (Al)]. As such h 
can be seen as a Hamiltonian in conventional QM, always with the same spectrum as H in 
PTQM. The associated eigenkets can also be mapped accordingly. Specifically, an eigenket 
\ipn) of H is mapped to an eigenket of h with the same eigenvalue, i.e., for = ri\ipn), 

H\4jn) = En\iJn) = S„|0n). (22) 

We next extend this mapping via t] to general time evolving states as well, i.e., 

\'^)=Vm, (23) 
where |\1') (|$)) represents the time evolving states in time-dependent PTQM (conventional 



QM). Using our general Schrodinger-like equation [Eq. (18)] for PTQM, one finds the fol- 



lowing equation of motion for a mapped state in conventional QM, 

ih^\^) = h\^) with h = h + — 'qrj'^ - {i)r]'^y . (24) 

Because 7777"^ is not Hermitian in general, h ^ h and the above mapped equation of motion 
is in general different from the standard Schrodinger equation in Eq. ^ for the mapped 
Hamiltonian h. 



Equation (24) indicates that among all possible choices for the mapping operator i], the 
one with 

^proper''7proper (^propor^proper y (25) 

is special because then, the mapped Hamiltonian operator h can afford to play a dual role 
as expected from conventional QM: It not only gives the spectrum of the mapped system, 
but also becomes the generator of time displacements for the mapped time evolving states. 
To emphasize such a peculiar type of mapping, we call this mapping a "proper" one. In 
this terminology, our time-dependent PTQM is equivalent to certain Hermitian problems in 
conventional QM under a proper mapping. Note however, it can be challenging to find this 
proper mapping explicitly, a strong indication that time- dependent PTQM is not a trivial 
extension of conventional QM. 
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To better understand the complexity of the proper mapping, we first discuss one possible 
procedure to find the proper mapping. Suppose an unknown ?7proper is a proper mapping 



satisfying Eq. ( 25 ) . Then another known mapping with 



rj = Urj 



proper 



is found to satisfy 



U 



rjTj 



nt 



U. 



(26) 



(27) 



Seeking the proper mapping %roper itself is now converted to solving the above first-order 
differential equation for U . This can be easier than directly searching for a Hermitian 



^proper 



proper 



. The suggested procedure hence goes like the following. First, we find a special 



solution to Eq. (19), denoted by rj. Then we solve Eq. (27) for U. Finally, -/^proper = U^il 



gives us a proper mapping. Because the right hand side of Eq. (27) is time-dependent, in 



general there is no analytical solution to U . In addition, there are no boundary conditions 



for U with respect to time [such as U\t=Q = 1], so the general solution to Eq. (27) has an 
arbitrary constant unitary factor. This factor corresponds to the unitary equivalence in the 
mapped Hermitian system, namely, the arbitrariness of choosing a basis to write down the 
Hamiltonian. 



For a proper mapping r/proper, the Schodinger-like equation in (18) takes the form 



^^di '^^ = ~ iHroperVoper) 1^)- 



(28) 



An equation of this form was also studied in Ref. [0]. However, we should stress that our 
perception and approach are distinct from Ref. [6]. In Ref. p], the author started from a 
Hermitian Hamiltonian h{t), and then applied a time-dependent similarity transformation 
to obtain H{t) = r]^^{t)h{t)r]{t) (Note that rj was called fl there) and the equation of 



motion Eq. (28). In such a treatment, a PT-symmetric system is derived from a Hermitian 



system (so everything is known from the Hermitian system). On the contrary, we treat a 
non-Hermitian "PT-symmetric Hamiltonian as a fundamental starting point, with a proper 
mapping between this system and a Hermitian system to be sought. 



Before ending this subsection, it is worth noting that Eq. (27) is also applicable to a 
Hermitian system in conventional QM. In that case, W = 1, and any unitary operator 



would be a solution to r] in Eq. (19) and the simplest solution r/propcr = 1 is a proper 
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mapping. Equation (27) then gives that an improper mapping r] shifts the Hamiltonian by 



rjT] 



(29) 



which is precisely the expected shift to the time-displacement generator in conventional QM 
under a time-dependent unitary transformation [as already seen in Eq. ([3])]. 



III. EXAMPLE: A COMPLEX HARMONIC OSCILLATOR 



In this section we wish to use a very simple example to illustrate time-dependent PTQM. 
Though in general it is challenging to identify how such time-dependent problems can be 
properly mapped onto a familiar problem in conventional QM, we have made efforts to 
construct examples where such a proper mapping can be found. 
Consider then a complex quadratic Hamiltonian 
1 



CHO 



Y 

X + 2iy- - y"— ) + (y + \y) {pq + 



+ Zp' 



(30) 



where X, Y, Z, and y are real, possibly time-dependent parameters, with X, Z, and y 
positive and ZX > Y^. It is straightforward to find that the eigenvalues of this complex 
Hamiltonian are all real and positive, which are given by 



E„ = (n + i) hVzx - r2 with 



n 



0,1,2,. 



(31) 



One also finds Hquo in (|30|) being self-adjoint with respect to the following metric operator 

(32) 



W 



exp 



hz 



According to our axiom proposed above, we then obtain the time-displacement generator A 
for this system. 



A 



y • d fy 



q^ + (Y + iy) {pq + qp) + Z'^p^ 



(33) 



To digest the above expression for A, we now attempt to find a proper mapping to identify 
an equivalent time-dependent problem in conventional QM. To that end we factorize W as 
W 



n 



proper ' 



?7proper- As sccu in Eq. (26), any other mapping rj can be expressed as a product 



of a unitary operator and the proper mapping, rj = f/?7proper- Using Eq. (32), one can easily 
find one Hermitian square- root of W, i.e.. 



?7 = exp 



'2hZ^ 



(34) 



To find a proper mapping, we then consider the following ansatz for the unitary factor, 



U = exp 



1 



(35) 



where k and v are real coefficients to be determined. 

The proper mapping operator ?7proper = U'^V maps Hquo onto a Hermitian Hamiltonian, 

^CHO = ''?proper-f^CHO'7proper 

= I [(X + + K^Z) f + (Y + kZ) {pq + qp) + Zp"] . (36) 

The mapped Schrodinger equation in conventional QM has the following time-displacement 
generator. 



'-CHO — ^CHO + |i^ 



rji] ^ ~ {r]rj ^)^ 



= hcno + + V. (37) 

Since ?7proper is a proper mapping, we must have /icho = ^cho- So k and d need to be chosen 
as constants in time. As a result, f/ is a constant unitary operator. For convenience we may 



choose K = and v = Q. In other words, the Hermitian mapping operator rj in Eq. (34) 
is actually a proper one. Finally, the time evolution of a complex harmonic oscillator is 
properly mapped to that of the so-called generalized harmonic oscillator in conventional 
QM, namely, 

/iGHO = \ [Xq^ + Y{pq + qp) + Zf] . (38) 

It should be stressed that this example is specially designed to demonstrate the consistency 
of time-dependent PTQM. In general it can be challenging to find the analytic form of a 
proper mapping. 



IV. BERRY PHASE IN TIME-DEPENDENT PTQM: HOW A CONTINUOUSLY 
TUNABLE FICTITIOUS MAGNETIC MONOPOLE EMERGES 

In this section we focus on adiabatic evolution and the associated Berry phase in a 
class of PT-symmetric two-level systems. When considering adiabatic manipulation, the 
Hamiltonian changes with time and so does the metric operator in general. As such, this kind 
of setup becomes a test bed to apply our time-dependent PTQM. Since Berry-phase problem 
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in conventional QM is well known, it will be also fruitful to make a comparison between 
predictions from time-dependent PTQM and the familiar results in conventional QM. In 
particular, we shall show how time-dependent PTQM may lead to surprising dynamical 
consequences and then explain them. 



A. Adiabatic evolution and the geometric phase in PTQM 

We consider a general time-dependent PTQM problem with the Schrodinger-like equation 



in Eq. (18), where is a time-dependent PT-symmetric Hamiltonian that determines the 
instantaneous spectrum of the system. The explicit time dependence of H is implemented 
via a set of time evolving parameters X = (Xi,X2,X3, • ■ ■ ) such that H = i/[X(t)]. The 
associated W can then be understood a.s W = W[K(t)]. 

To investigate the adiabatic evolution, let us first take a time derivative on the eigenvalue 
equation ([s]), 

H\^n) + H\ijn) = Enl^n) + ^n|^n). (39) 

Multiplying the above equation by from the left and then applying the orthonormal 

condition in Eq. (|6]), we arrive at 

{ij^\WH\iJn) = 5mnEn + - ^„) | | . (40) 

Assuming that there is no degeneracy in our Hamiltonian, the above equation leads to 

»J^^'W.„) = rn^n, (41) 

and 



En={^n\WH\^n). (42) 

Next, we expand in terms of the eigenstates of H, 

|vl/) = 5^c„e^""|V^„), (43) 

n 

where a„ is a dynamical phase defined as a„ (t) = -|/ drE„[X(r)]. Taking the time 
derivative of this expansion we obtain 

1^) = 5Z - ^^"^") ^""1^") + J2^ne"'"\i^n). (44) 

n ^ ^ n 
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Plugging the above expression into the Shrodinger-hke equation in (18) and applying the 
eigenvalue equation in ([s]), we get 

n n n 

Further multiplying the above equation by {ipmlW from the left and then using the relations 



in Eqs. (41 ) and (42 ), one finds the time dependence of the expansion coefficients Cm, namely, 



Cm = -Cm { (ipmlWl^Jm) + ^(^m|W^|^m) 



n^m \ / 

As in conventional QM, we now make an adiabatic approximation by ignoring the population 
transition terms (i.e., with n 7^ m) in the above equation. This then brings us to 

c™(t)~c^(0)e^^'"W, (47) 

where 7^ is identified as the geometric phase and it is given by 

iL = i {{i^m\W\i;m) + l{i^m\W\^m)) • (48) 

In terms of time evolving system parameters X, the geometric phase is then given by 

7^ = iy"dX- [{^m\WV\^m) + lii^mliVWMm)] , (49) 

where V = The geometric phase associated with a closed adiabatic path yields a Berry 



phase in our time-dependent PTQM. Interestingly, the first term in Eq. (48) or Eq. (49) 
represents an integral over a Berry connection (as expected, the connection as an inner 
product carries a metric W). On top of that, there is a second term, which originates from 
the metric's dependence on system parameters. Only the combination of the two terms 
results in a real phase 7^. 

B. Berry phase in PT-symmetric two-level systems 

Exploiting a full parametrization of 2 x 2 PT-symmetric Hamiltonian matrices we 
have the most general Hamiltonian of a 2 x 2 PT-symmetric system, 

H2X2 = el2x2 + (a n"" + ib sin Sn'^ + ib cos 5 n"^) • cr, (50) 

12 



where all six parameters, e, a, b, 6, ip, and S are real, cr are Pauli matrices and three mutually 
orthogonal unit vectors are defined as 

n^' = (sin 6' cos v?, sin 6' sin cos 6'), 

= (cos 6 cos (f, cos 6 sin (p, — sin 6) , 
n*^ = (— sin(/9, cos 0). (51) 

Obviously, i?2x2 is Hermitian iff 6 = 0. 

For each angular parameter, 6, (f, or S , H2X2 is a periodic function with the same period 
27r. Since a — )■ —a is equivalent to — t- ± tt and 5 — ?■ 27r — 5, we may work with a > 
without loss of generality. This choice will help us to simplify the expressions of eigenstates. 

The eigenvalues of -^2x2 are 

E±=e± Va2 - 62. (52) 

They are real when a'^ > b"^. For our later calculations, we limit our discussions in the regime 
of > 6^ where -^2x2 is always diagonalizable with two non-degenerate real eigenvalues E± 
(PT-symmetry is not spontaneously broken). The corresponding eigenstates are 

(cos f + I3e-'^ sin |) e"'"^ 
-/^e^^"^ cos f + sin I 

- {l3e'^ cos f + sin I) e" 
cos f - f3&^ sin f 

(53) 




where B = — K „ . In the Hermitian limit, = 0, the above eigenstates reduces to the 



conventional choice for the Hermitian 2x2 Hamiltonian [see Eq. (75)] 



The metric operator W is found to take the following general form, 

W = /x[al2x2+ {un^ + bcos6n^ -bsm6n^) - a] , (54) 



where fx is positive and u is sufficiently small, \u\ < ^ — h^. The eigenstates in Eq. (53) 
are orthogonal with respect to this metric operator, 

{%l)+\W\tl).) = 0= {i)4W\^^+). (55) 

The norms of the eigenstates are 



A4 = = 2/x |n±|^ y^^TT^^^^^^i^. (56) 

va — + a 
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Note that A/± are always positive. That is, the eigenstates are always normalizable in 
the above ranges of parameters. To normalize the eigenstates, one may either adjust the 
coefficients n± in \4'±), or tune parameters /i and u in W. These two ways are equivalent. 
This means that the undetermined parameters fi and u can always be absorbed by the 
normalizations of the eigenstates. We choose fi = 1/a and i/ = for our considerations 
below. By fixing /i and u in this manner, the parameters in W are the same ones in H, i.e., 
W = W[K{t)]. With these simplifications, W assumes the following form 

W = l2x2 + ^ (cos 5 - sin 6 n'^) ■ cr. (57) 

In the Hermitian limit, b ^ 0, W reduces to the identity matrix which is the conventional 



choice. To normalize the eigenstates using the metric in Eq. (57), the coefficients in Eq. (53) 
can be chosen as 

"^ = ^ 2(a^-6^) ' 

where we have explicitly chosen one realization of the rather arbitrary overall phase factor 
to ensure that the eigenstates are periodic function of 9 with a period 2tt instead of 47r. 



Taking note of both expressions of -^2x2 in Eq- (50) and W in Eq. (57), we shall consider 
adiabatic processes in which the three angular parameters (6', (p, 5) change with time and an 
adiabatic cycle is formed in the end. 

Using general results in the previous subsection, lengthy calculations eventually lead to 
the following surprisingly simple expressions for two geometric phases 7^ associated with 
two adiabatic eigenstates 

1± = j d(/^ + Flde + Fi d5) , (59) 

with 

F^ = - ( l±—^=cos9] , 
F<^ - ^ 

Note that the expression of F^ depends on the phase factor in the normalization constants 



n±. For example, if we choose the factor to be e'2 in n±, then we would have F^ 



1 

± - 2- 
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However, the difference is not detectable for 9 varying from to or to 271. We also emphasize 



that the integral in Eq. (59) is to be carried out along a line defined by the adiabatic path 
in the parameter space of {6, ip, 5). 

One may notice that the results of and here differ from our earlier results obtained 
in Ref. [1] [Eqs. (23) and (24)]. This is due to a combination of two different choices we 
make here as compared with Ref. The first different choice is in the phase factors of 



the eigenstates. Our choice here in Eq. (53) guarantees that both eigenstates are periodic 
function of 9 with a period 27r. By contrast, in Ref. the parameter 9 was always treated 
as one of the spherical coordinates varying from to tt, so the eigenstates in Ref. |1] [Eq. (19)] 
need not to be, and were not chosen to be, periodic function of 9 with period 27r. The second 
different choice is in the time-displacement generator A. Using the current language, the 



choice made in Ref. [1] accepts H ^ H, which is no longer advocated here. 

Note that = = 0, which implies that for an arbitrary but fixed y?, there exists 
full integrals G±{9,6) such that 

dG±{9, 6) = Fld9 + Fl dS. (61) 

Indeed, we can easily obtain 

Because G±{9,6) is not a periodic function of either 9 or 6, a path along which both the 
Hamiltonian and the metric have returned to their initial forms (hence regarded as a closed 
path) may yield a nonzero Berry phase. For example, after 6 varies from to 27r with both 
9 and fixed, the Berry phase is found to be 

7| = G±{9, 2vr) - G±{9, 0) = ± (l - ^7^^) vr. (63) 

Similar considerations apply to the parameter pair (yj, 5) because -q^ = = 0. We 
will not repeat the analogous calculations here. 



Just like in conventional QM, cases with time- varying parameter pair {9, (p) turn out to 

dFt , dFl 
dip 



be more involving and interesting. In that case, -^^ 7^ -jr^ and hence there are no simple 



full integral solutions for the geometric phase. 
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C. From Berry phase results to a continuously tunable magnetic monopole 



For Berry phase generation via the adiabatic parameter pair {9, (p), we adopt a language 
similar to our previous work To proceed we assume {6, (p) are the two angles in a standard 
spherical coordinate system, (r, 6', yj). We still denote 7^ as the Berry phase generated after 
the system has adiabatically moved along a closed path on a sphere, with 

li = j (^±d(^ + Fide) . (64) 

Next we define the following vector potential field, 

A± = - f^n^ + ^n-^ , (65) 



e \ r rsin^^ 

where e has the dimension of an electric charge. It is straightforward to verify that 

Fide + Fidyp = ■ dr. (66) 



One therefore has 

/ A± ■ dr. (67) 



hi 

Let us then consider the surface enclosed by a closed adiabatic path on the sphere. 
If it does not include the north pole, then one finds that the vector potential field A is 
well behaved everywhere on this surface. As a result, one may apply Stokes' theorem and 



reexpress Eq. (67) with the surface integral 

7f = ^ V X A± ■ dS. (68) 

Further using = and the explicit form of the curl operator in a spherical coordinate 
system, this surface integral reduces to 

where Vt is just the solid angle covered by the adiabatic path on the sphere. 

On the other hand, if the adiabatic path on the sphere does enclose the north pole, 
then the vector potential field A diverges at ^ = and hence it is illegitimate to apply 
Stokes' theorem directly. One simple remedy is to divide the adiabatic path into two parts: 
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the first part is a deformed patli circumventing tlie nortli pole and tlie second part is an 
infinitesimal circle surrounding the north pole. In this manner, one may apply Stokes' 
theorem again to the first part and directly calculate the line integral for the second part. 
Without loss of generality we always assume counter-clockwise adiabatic paths. Then the 



first part contributes to the line integral in Eq. (67) by again ^^^^^==^fl, and the second 
part contributes 2TTF^\g=Q = ^1 ± yj^pj tt. Summing up the two contributions, we find 

if a closed adiabatic path encloses the north pole. 

In conventional QM the Berry phase of a two-level system in a rotating field can be 
interpreted as the flux of a fictitious magnetic field. We may do the same here. In partic- 
ular, the above expressions of 7^ can be related to a solid angle Q, plus a constant term 
depending upon if the adiabatic path encloses the north pole. This suggests that 7^ may 
be connected with the flux of a magnetic monopole plus a singular field along the ^ = 



direction. Specifically, by combining Eqs. (69) and (70), one obtains the following unified 
Berry phase expression. 



7i = ^yyB±.ds, (71) 

where the fictitious magnetic field is given by 

„ Tih f ^ , a \ ^, . ^ h a r 

with = n''|e=o- This virtual magnetic field describes a monopole with charge 
^J^i,i Qrri^ at the origin and a string component along the 6* = (2;-) direction, where 
= ^ is the conventional magnetic monopole charge [B]. In sharp contrast to the 
well-known Berry phase problem of a two-level system in conventional QM, here the ficti- 
tious magnetic monopole is no longer quantized: its charge is continuously tunable between 
(—00, —Qm\ and [Qm-, 00). Echoing with this, the phase contribution from the string com- 
ponent is no longer an integer multiple of 27r and hence becomes an observable quantity. 
This is again drastically different from the trivial Dirac string in conventional QM, which 
can only produce an unobservable phase. 
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D. Proper mapping to Hermitian two-level systems 



The Berry phase results in the previous subsection are surprising and call for more analysis 
and theoretical insights. Consider then a proper mapping from the above PT-symmetric 
two-level problem to a Hermitian one in conventional QM. We proceed by first choosing a 



Hermitian similarity transformation operator 77, the square-root of W. For W in Eq. (57) 
it has two distinct square-roots (on top of the overall ± signs) [7], 

1 / y'^^'') — 6sin6'cos5 6(cos6'cos5 + isin5)e~'''' , 
V^^ = I (73) 

V 2ax(±) \ 6(cos 6* cos (5 - i sin (5)e"^ x^"^^ + 6 sin 6^ cos 6 



with x^^^ = a±-\/a^ — b'^. To avoid confusion, we always use subscripts ± to label eigenvalues 
E± and use superscripts '•^^ to label the different choices in the mapping r/^^^ 

Under the similarity transformation t]*-^-*, we have two surprisingly simple Hermitian 
Hamiltonians 



-1 



= el2x2 ± y/a^^-l^ if ■ cr 

e o\ . [ cose sin^e"^^ , , , 

± - 62 . (74) 

e / I sin^e"^ -cos^ 

Note that h^'^^ ^ h^~'^ when 6 is shifted by tt, namely, /i*-'*'^ and h^~^ are actually the same 



Hamiltonian with a shifted parameter 9. Since the ± signs in front of ^0? — b'^ in Eq. (74) 
swap two eigenvalues, all the eigenstates of h'^'^^ are also eigenstates of /i^^^ with swapped 
eigenvalues. Namely, the normalized eigenstates of /i'-^-' are 



cos I e 
sin ^ 



f:\ r \ .9 / — sin^e , 

^(+)) = 10^-)) = e-'5 \ , (75) 



cos ^ 

where we have made an overall phase convention as in the PT-symmetric case. The advan- 
tage for this choice is that in the limit of 6 0, we have 

Iimi72x2 = and lim l^-i) = 10^+^). (76) 

6— s>0 b^O 
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Neither 77*^+^ nor i]^ Ms a proper mapping as defined in Eq. (25). That is, the mapped 
equation of motion will not be governed by h^'^\ To find a proper mapping, we need to 



solve Eq. (27) for the unitary matrix U which links a proper mapping proper to the current 
known (improper) mapping t]^"^^ by i]^^-* = f/r/proper- For an arbitrary path with varying 9, 
(fi, and S, analytically obtaining such U to establish a proper mapping is impossible. Below 
we consider paths that involves changes in just one of the three angular variables. 



1. Varying 6 only 



Plugging the improper mapping in Eq. (73) into Eq. (27) and keeping in mind that only 



5 is a function of time, we get a differential equation for t/ as a function of 6, 



a6 where C^"^ - 2 ^ ^^=J • (77) 

The ranges of C^^'' are C^'^^ ^ and ^-^^ > 1, where the equal signs are only taken at the 



Hermitian limit, b ^ 0. Equation (77) is easy to solve since the factor in front of U{6) on 



the right-hand side is independent of 6. The solutions are 

U{6) = exp(iC(5n'' ■ (t)Uo 

cos{(6) + i cos 9 sin(C5) i sin 9 sin(^5)e~"^ 
i sin 9 sin(^(5)e"^ cos{(6) — i cos 9 sin(^5) 



Uo, (78) 



where Uq is an arbitrary constant 2x2 unitary matrix. For a shorter notation, sometimes 
we suppress the superscripts in (^"^^ 

Remarkably, though the PT-symmetric Hamiltonian H2X2 is a periodic function of 6 
with a period 27r, the found proper mapping, i.e., U'^{5)r]^'^\ is not a periodic function of 
6 in general. As a consequence, when the parameter 6 varies from to 2-7?, the mapped 
system does not return to itself in general. This explicitly demonstrates that mapping one 
time-dependent system to another must be carried out carefully. 

Another interesting observation can be made. The above matrix structure of exp{i(Sif -cr) 



is much analogous to that of /i(^) in Eq. (pi. The consequence is that the two matrices 



commute with each other. As such, when only S varies with time, the properly mapped 
Hermitian Hamiltonian is simply given by Uoh^^^Ul where Uq is a constant and arbitrary 
unitary matrix. This being the case, the parameter 6 completely drops out from the mapped 



19 



Hamiltonian. However, the parameter S still appears in the mapped states. For simplicity, 
let us take Uq = ^2x2 and only consider the mapping to h^'^\ In this case, the proper 
mapping rjptoper = [U^~^\6)y r]^'^^ maps the eigenstates in Eq. (53) to 

vitU^±)=e'^-''^'^V±), (79) 

where 10^"^) are the conventional choice of the eigenstates for the Hermitian Hamiltonian 
h^^^ in Eq. (75). To understand the extra phase factor, e^'"''^''^, let us recall the calculation 
of the Berry phase in a "PT-symmetric system. In the adiabatic limit, if we start with 



an eigenstate of -^2x2, the wavefunction [a solution of the Schrodinger-like equation (18) 
remains as an instantaneous eigenstate. The phase of the wavefunction has two parts, a 
dynamical phase and geometric one. For example, 

^e'^±+'"±|^±), (80) 



where 7^ is the geometric phase given in Eq. (59) with d6 = dip = 0, a± is the dynamical 



phases, and \ip±) are eigenstates given in Eq. (53). Using the proper mapping, we obtained 
the mapped time-evolving wave function at the end of the adiabatic process, i.e., 

\^±) = vitU'^i) ^ e'^±+-±^^^'"'V^V (81) 

Since the mapped Hamiltonian h^~^^ is time- independent, the phase of the mapped wave- 
function must be the dynamical phase and nothing else. Therefore, we get 

li = ±C(^^5 = ±Ul- -^4=^) S, (82) 



2 V Va^ - 



this is precisely the same geometric phase one would directly get from Eq. (59) if only 5 
varies with time. The rather strange Berry phase result in the simple PTQM example here 
is thus well understood. In general, when 6 or (p also changes with time, the proper mapping 
is much more complicated and cannot be reduced to a simple phase factor. 



2. Varying 9 only 



In this case, Eq. (|27|) reduces to 

sin 6* sin 5 (— cos 6* sin 5 + i cos 5)e~"^ 

cos0sin(5 — icos(5)e^''' — sin ^ sin 5 



dU{9) 
dO 



i( cos 6 



u{e) 



-iC cos(5e-^2'^3e-'2''2e'2''V2e-^2'"3gi2'^2gif^3^(^^)_ 



(83) 
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This equation is not easy to solve because the prefactor of U{9) on the right-hand side 
depends on 6. Noting the factorization of the prefactor, we consider a differential equation 
for 

v{e) = S^'S^^u^e), (84) 

which is independent of 6, i.e., 



^ = i [-Ccos^e'i-V^e-'i'^^ + \a, ) V{e). 



(85) 



This differential equation can be integrated directly. In the end, the solutions to the original 



differential equation in Eq. (83) are found to be 



u{e) 



e-if '^3e"'5'^2 exp 



i (-Ccos(5e^5-3^2e-'^-3 + i^^) 6 



(86) 



From this expression for U{9), the obtained proper mapping W{9)ri^'^'' is not a periodic 
function of 6 in general. As to the properly mapped Hermitian Hamiltonian, we parameterize 
it by 

cos[e(e)] sin[e(^)] e-'*^'') 
sin[e(^)] e'*(^) -cos[e(0)] 
For simplicity, we choose Uq = l2x2- Then the angular parameters for h{9) are found to 
satisfy 



(87) 



cos[0(^)] = cos eVl-4C(l-C) cos2(5 
/l-C(l + e 



(88) 



1 - C (1 + e^i-^) ■ 

Interestingly, this mapped Hermitian Hamiltonian does not depend on ip. This is accidental 
due to the arbitrariness in choosing Uq. For example, for an alternative choice of the constant 



unitary matrix given by f/^ = e^f'^3^ e 

remains the same as in Eq. (\88h but $ is now 

determined by 

'l - C(l + C-2i5) 



l-C(l + e 



(89) 



To solve for Q and $ from Eq. (88) or Eq. (^89^, the multi- valued inverse trigonometric 



functions may be used. In doing so one should be careful in choosing the right branch of 
the inverse functions. 
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3. Varying if only 



In this third case, Eq. (27) gives 

cos^ 9 cos^ 6 + sin^ 6 



dip 



sin 9 cos 5(cos 9 cos 5 + i sin 6)e 



sin6'cos5(cos6'cos5 + isin(5)e^'^ — 

— cos 9 — i sin 9 sin 5 
i sin 9 sin 5 cos 9 



iC e-^2'"3e-'2'"2e'2''^ 



cos^ 9 cos^ (5 — sin^ 5 



(90) 



Use a technique similar to the previous case, we obtain the solution for U{ip) 

— cos 9 



U{(fi) = e-'^"^ exp <^ i 



isin6'sin5 \ .g 1 

i sin ^ sin 5 cos^^ ' ^ 



(91) 



For generic values of 9 and 5, the mapped Hermitian Hamiltonian is too complicated to be 
written down here. Instead, we only give the explicit forms for the special cases of 5 = 
and 5 = 71/2 and choose to numerically plot a generic mapping in Fig. [l| and Fig. |2| For 
simplicity, we choose Uq = 12x2 in the mapping. 



e / \ sin^ei(i-2C)^ 



(92) 



6 = 7t/2. In this case, the mapped Hamiltonian is very simple: 

eO\ ^/ cos^ sin^e-'^i-^^)^ 

— cos 9 

It is clear how the proper mapping changes the periodicity of the system. For a non- 
integral 2^, the mapped Hamiltonian does not come back to its original value when ip 
changes from to 27r. Since C^^^ < and C^^^ > 1, the factor |1 — 2C| is always greater 
or equal to one. This means that the period of mapped Hamiltonian, 2tt / \1 — 2(\ is 
always less than or equal to that of the PT-symmetric system. 

6 = 0. We may again use G and $ to parameterize the mapped Hermitian Hamiltonian, 



as in Eq. (|87|). In this case, the angular parameters in the mapped Hamiltonian are 

2C 



cos[9(v9)] 



cos( 



1- ^sin2e[l-cos(e<^)] 



sin[e((^)]e-^*(^) = sin^jl - ^ (1 -2Ccos2^) [1 -cos(^(^)] - ^sin(e(^)| 



(93) 
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where ^ = y^l + 4:( {( — 1) cos^ 6. The same key feature of the mapped Hamiltonian 
appears: Though it remains periodic functions of (f, it has a period different from 27r. 
Since ^ > I, the mapped Hamiltonian always has a shorter period, 27r/^. 

• Generic S. The expressions of the mapped Hamiltonian is too complicated to be 
included here. Therefore, we solve the angular parameters B and $ numerically for 
a special case with (^'^^ = —0.95, 9 = 7r/3, and S = 107r/9 and plot the results in 
two figures. In Fig. [T| we plot $ as a function of (p. It clearly shows that $ increases 
faster than (p. In particular, $ changes from to 27r as (p changes from to l.OTvr. In 
Fig. [2| we plot two mapped paths on the same unit sphere. The {9, Lp) path is plotted 
with the dashed red line and the (O, $) path is plotted with the solid blue line. The 
{9, (fi) path is horizontal because 9 is fixed; whereas the (O, $) path is not horizonal 
in general as both and $ change with ip. As the (G, $) path completes a circle on 
the sphere, we see that the parameter (f only changes from to I.OTtt, thus forming 
an open path instead. The proper mapping between an adiabatic process in PTQM 
and an adiabatic process in conventional QM is thus subtle. 

FIG. 1: Angular parameter $ in the mapped Hamiltonian as a function of the angular parameter 
(p in the "PT-symmetric Hamiltonian H2x2- Here we choose parameter ^(+) = -0.95 and fix the 
two other angular parameters at 9 = tt/S and 5 = lOvr/O. This plot shows that $ grows faster 
than ip and ^{ip ~ l.OTvr) = 2tt. 

4 

3 

1 


0.0 0.5 1.0 1.5 :.o 
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FIG. 2: The adiabatic path in PTQM and the mapped adiabatic path in conventional QM, plotted 
on the same unit sphere. The red dashed line is an open path in (9, ip) parameter space of the 
"PT-symmetric Hamiltonian H2x2 with 9 = ir/S and (p varies from to l.OTvr. The blue solid line 
is the mapped closed path in (G, $) parameter space of the mapped Hamiltonian. In the plot we 
choose parameters C^^^ = —0.95 and 6 = lOvr/Q. 




Summarizing the three cases analyzed above, it is clear that the properly mapped Hamil- 
tonian in general is not a period-27r function of the angular parameters 6, ip, or 6. From 
this viewpoint, our previous Berry phase results can be better understood. The Berry phase 
for H2 x2 in our PTQM is the geometric phase associated with a closed path in the [6, ip, 6) 
parameter space, which is in general an open path in the properly mapped (0, $) space. 
Hence, in the language of a properly mapped Hermitian two-level Hamiltonian [e.g., h{6) 



in Eq. (87)], our Berry phase results in PTQM problems can at most be translated into a 
geometric phase associated with an open path in the (G, $) parameter space. Equivalently, 
a closed path in conventional QM is mapped to an open path in PTQM. Therefore, there is 
nothing inconsistent between our axiom for PTQM and conventional QM. It is also clearer 
by now, that despite the existence of a proper mapping from time-dependent PTQM to 
conventional QM, time-dependent PTQM still has its own rich features. 
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V. CONCLUSIONS 



In this work we have advocated a general framework for time-dependent PTQM and dis- 
cussed in depth its relation with time-dependent problems in conventional QM. Our study 
is complementary to current studies of spectral problems in PTQM. Considering recent ex- 
perimental interests in PTQM pHTS]. two types of experiments on time-dependent PTQM 
may be called for. The first type is to simulate or synthesize the Schrodinger-like equation 
proposed in this work and then examine the time evolution, which is unitary if viewed with 
an appropriate time-dependent metric operator. In particular, the confirmation of our Berry 
phase results, which may be interpreted as the fiux of a continuously tunable fictitious mag- 
netic monopole plus a singular string component, would be of great interest. The second 
type of experiments is to verify our axiom for time-dependent PTQM, by (i) directly work- 
ing on certain non-Hermitian quantum systems with a time-dependent instantaneous real 
spectrum and a time-dependent metric, and (ii) checking if our proposed Schrodinger-like 
equation describes the actual quantum dynamics. Such experiments would constitute a truly 
fundamental test of time-dependent PTQM and they are not expected to come very soon. 
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Appendix A: The Uniqueness of the T^T-symmetric and "PT-anti-symmetric Parti- 
tion 

If an operator H is PT-symmetric with respect to W = rj^r], then rjHr]'^ must be 
Hermitian: 



-1 



(Al) 



Similarly, for a PT-anti-symmetric operator A, rjArj ^ is anti-Hermitian: 



WA = -A^W, r]Ar] 



-1 



(A2) 



Now, suppose that A has two different partitions. 



A = H + A = H' + A'. 



(A3) 
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Let us sandwich the above equation between t] and ri~^, obtaining 

■qHr]'-^ + r]Ari'^ = riH'ri'^ + rjArj'^. (A4) 

Note that both rjHri'^ and rjH'r]'^ are Hermitian and both rjAr]'^ and r]A'r]~^ are anti- 
Hermitian. Because the way of partitioning a matrix into a Hermitian part and an anti- 
Hermitian part is unique, we must have 

rjArj''^ = r]Ar]'^. (A5) 

This means 

H = H', 

A = A'. (A6) 

That is, the way of partitioning A into a PT-symmetric part and a PT-anti-symmetric part 
is unique for a given W. 
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